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A black hole immersed in a rotating Universe, described by the Gimon-Hashimoto solution, is 
tested on stability against scalar field perturbations. Unlike the previous studies on perturbations 
of this solution, which dealt only with the limit of slow Universe rotation j, we managed to separate 
variables in the perturbation equation for the general case of arbitrary rotation. This leads to 
qualitatively different dynamics of perturbations, because the exact effective potential does not 
allow for Schwarzschild-like asymptotic of the wave function in the form of purely outgoing waves. 
The Dirichlet boundary conditions are allowed instead, which result in a totally different spectrum 
of asymptotically Godel black holes: the spectrum of quasinormal frequencies is similar to the one 
of asymptotically anti-de Sitter black holes. At large and intermediate overtones N, the spectrum 
is equidistant in N. In the limit of small black holes, quasinormal modes (QNMs) approach the 
normal modes of the empty Godel space-time. There is no evidence of instability in the found 
frequencies, which supports the idea that the existence of closed time-like curves (CTCs) and the 
onset of instability correlate (if at all) not in a straightforward way. 



PACS numbers: 04.70.Bw,04.50.-h,04.30.Nk 



I. INTRODUCTION 

Various structures in the universe starting from com- 
pact objects, such as stars and planets, and including 
constellations, galaxies and clusters of galaxies, rotate. 
This made physicists suppose that the Universe itself can 
be rotating. The first model for such rotating universe 
was suggested by K. Godel yet in 1949, who found an ex- 
act solution for the 3 + 1 dimensional rotating Universe 
[l| in General Relativity. The solution is homogeneous, 
has rotational symmetry, and allows for the definition of 
the direction of positive time consistently in the whole 
space-time, and, what has been in the focus of recent 
research. It also allows for closed time-like curves. 

In the Einstein gravity a static black hole cannot be 
embedded in the vacuum Godel Universe. However such 
an embedding is possible if extra fields are added. The 
asymptotically Godel, 5-dimensional, nonrotating black 
hole solution was obtained by Gimon and Hashimoto Q 
for the low-energy supergravity limit, i. e. in the pres- 
ence of an extra gauge field. This solution attracted con- 
siderable interest due to its utility when studying exact 
solutions in the form of the so-called pp-waves in string 
theory ||. 

Here we are interested in perturbations of asymptot- 
ically Godel space-times from a different point of view. 
First, the Godel Universe has a peculiar feature: at its 
asymptotic, far region it allows for the closed time-like 
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curves (CTCs) which opens possibility of existence of the 
time machine. Thus, existence of compact objects (for in- 
stance, black holes) in such space-times would look most 
unreal and one could expect that there should be some 
physical laws which forbid such geometries. Indeed, there 
are a number of topological theorems which shows that 
nature should favor "normal" causality of the world. Yet, 
these are rather geometrical arguments, than a rigorous 
physical reason. 

An interesting observation, which might be an exam- 
ple of such a physical reason of the nonexistence objects 
allowing for CTCs, was made in where it was numeri- 
cally shown that the infinite cylinder (with an event hori- 
zon) in General Relativity has infinitely growing modes in 
its characteristic spectrum. Thus, the cylindrical space- 
time, allowing for CTCs, is also gravitationally unstable. 
From this observation, one could expect that there is 
a deeper correlation between the onset of the instabil- 
ity and existence of CTCs. A straightforward way is to 
check other compact objects with CTCs as to their sta- 
bility. This can be done by investigation of the proper os- 
cillation frequencies (called quasinormal modes [5j-|7| for 
black holes) of the compact object under consideration. 
These frequencies dominate at intermediate late time in 
the black hole's response to the external perturbation. 

Quasinormal modes of the Gimon-Hashimoto black 
holes were first studied in [8], though the wave equa- 
tion in [8j was obtained for the slow rotation limit. In 
this limit, the obtained effective potential looks like the 
one for the 5-dimensional Schwarzschild black hole with 
a shift in frequency, which depends on the Universe ro- 
tation parameter j. In this approximation properties of 
the asymptotics of the wave equation and, consequently, 
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boundary conditions, are qualitatively the same as for 
the Schwarzschild black hole. Therefore, quasinormal 
spectrum is just shifted Schwarzschild spectrum H . Fol- 
lowing the approximation used in [8j], there appeared a 
number of consequent works on quasinormal modes and 
Hawking radiation of more general solutions and other 
fields j|. 

Thus, we have a few motivations to study QN spec- 
trum of the asymptotically Godel black holes: First is 
the revealing of possible correlation between the gravi- 
tational instability and existence of CTCs, second, quest 
for a qualitative insight of main features of the quasinor- 
mal spectrum in the rotating Universe, and third, possi- 
ble utility of QNMs of the exact sting theory background 
given by the Gimon-Hashimoto solution. Having in mind 
the above motivations, in the present work we managed 
to separate angular variables without any approximation 
through using different coordinate system. The exact 
wave equation which we obtained revealed qualitatively 
new feature, which is connected with different asymp- 
totics of the wave equation in the region far from a black 
hole. We have found that no Schwarzschild-like boundary 
conditions, that is, purely outgoing waves, are allowed at 
infinity. Instead Dirichlet boundary conditions should 
be imposed, in a similar fashion with the asymptotically 
anti-de Sitter black holes, whose QNMs were extensively 
studied in the context of gauge/gravity duality (see [6| 
and 0] for recent reviews of QNMs). 

The paper is organized as follows. In Sec. II we men- 
tion some basic properties of the Gimon-Hashimoto solu- 
tion. Sec III describes separation of variables in the wave 
equation for a scalar field. Sec. IV is devoted to bound- 
ary conditions for the quasinormal modes, while in Sec. 
V we briefly relate the numerical methods which were 
used for calculations of the quasinormal modes. Sec. VI 
discuss particular case of empty Godel space-time (with- 
out a black hole) . Sec VII and VIII are devoted to spec- 
trum of Godel black holes, which includes the case of 
the quasiextremal Universe rotation for which some an- 
alytical results were obtained. Finally, we discuss the 
obtained results and their meaning for stability of the 
black holes and its possible correlation with the existence 
of closed time-like curves. 



II. 5-DIMENSIONAL ASYMPTOTICALLY 
GODEL BLACK HOLE: THE 
GIMON-HASHIMOTO SOLUTION 

The bosonic fields of the minimal (4+1)- supergravity 
theory consist of the metric and the one-form gauge field, 
which are governed by the following equations of motion 
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(i) 

(2) 



In the Euler coordinates (t,r,6,ip,(j>), the solution of 
the equations of motion JT][5]), describing the Godel uni- 
verse, has the form 

ds 2 = -(dt+jr 2 a 3 L ) 2 + dr 2 

2 

+ V -^(d6 2 + dip 2 + dcf 2 + 2cos0##), (3) 

where a\ — d<p + cosOdijj. The parameter j defines 
the scale of the Godel background. At j = we 
have the Minkowski space-time. The solution for the 
Schwarzschild black hole in the Godel universe is given 
by 1 

ds 2 = -f(r)dt 2 - g(r)ra 3 L dt ~ h(r)r 2 (a 3 L ) 2 + k(r)dr 2 
r 2 

+ — {d9 2 + d^ 2 + d4> 2 + 2 cos 9d^d4>) , (4) 



where 



g(r) = 2jr, h(r) = f(r 2 + 2M), 
f(r) = l-— , k(r)= 1-— + — V- 



(5) 



The radius of the event horizon is also corrected by pa- 
rameter j, 



r = y/2M(l-8j 2 M). 



(6) 



An essential feature of the above Gimon-Hashimoto so- 
lutions is possibility of the closed time-like curves (CTCs) 
in a region far from the black hole, 



r > 



y/l ~ 8Mf 
2j 



(7) 



III. SEPARATION OF THE VARIABLES 

In order to derive the wave equation we only need the 
relation 



g 2 (r) + f(r)(l-4h(r)) 



1 



jfe(r) 



which implies 



detg^ 



r 6 sin 2 1 
64"" 



The nonvanishing components of the inverse metric 
have the following form 

g" = -(l-4h(r))k(r), ff pr = fc- 1 (r), g ee = 4/r 2 , 
^ _ 4 ^ _ 4(cot 2 fl + fc(r)/(r)) 

r 2 sin 9 r z 
M - 2 9{r)k{r) m _ 4cosfl 
r r 2 - sin a 
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The scalar field perturbations in a curved background 
are governed by the Klein- Gordon equation 

= {fy/^gdvti) = m 2 *- (8) 

Since the background metric has the Killing vectors dt, 
d^, dij,, we choose the ansatz for the wave function as 



_ „—iujt + imp + im 



(9) 

Substituting ([9]) into ([H} and separating the variables, 
one can find that the angular part of the function satisfies 
the equation 

1 d . . d 2mn cos 9 — to 2 — n 2 ,\ 

■smO— H h A Y(6) = 0, 



sin 9 d9 



d6 



(fO) 

where A is the separation constant with the eigenvalues 

\ = 1(1+1), £ = max(|m|,|n|)+i, i = 0,1,2.... 

The equation for the radial part takes the wave-like 
form 

/ ^ + O(»-*))«(r*)=0, (11) 

where r* is the tortoise coordinate, which is defined as 

dr* = k(r)dr. (12) 

The effective potential can be written in terms of the 
coordinate r as follows 



Q(r) = (1 - 4h(r)) to 



2mg(r) 



r(l - Ah{r)) 



(13) 



1 /4A 



k(r) \r 2 



I6m 2 h(r) 3fc(r) — 6rfc'(r) 



(1 - 4/i(r))r 2 



4r 2 fc 2 (r 



Now we are in position to perform numerical analy- 
sis of the quasinormal spectrum for the obtained wave 
equation. 

IV. BOUNDARY CONDITIONS 

At the classical level, a wave cannot be emitted from 
the black hole horizon, so that the standard boundary 
condition for a great variety of black holes is requirement 
of the purely ingoing wave at the event horizon, 



Roc e-^V 1 - 4 ^''*) 



r* -» -oo, (14) 



where 



2mg(r ) 



r - 4r /i(r ) 
Taking account of (fc _1 (ro))' = 2/ro, we find that 



Roc(r-r )-^ ro ^ 1 - 4h(ro)/2 , r -> r„. (15) 



At spatial infinity we have 

i?(r -> oo) = C+*+(r) +C_*_(r), (16) 

where 

M± , A 2± , A 3 ± 



*±(r) = e ±Jwr r Q± 1 



y 3 



with 



- ± I 2m + QujjM - 32wj 3 M 2 



io 2 - (i 2 



_l_-2 

Since the exponents e JUJT have purely real index, 
they do not describe ingoing or outgoing waves. There- 
fore, we are unable to impose usual quasinormal bound- 
ary conditions. However, we can use the analogy with 
AdS backgrounds and require Dirichlet boundary condi- 
tions at spatial infinity. This implies that C + = for 
Re(ju) > or C_ = for Re(juj) < 0. 

The nontrivial behavior of the functions ^S± is observed 
when Re(uj) — 0. In this case both exponents have os- 
cillatory behavior at spatial infinity. Thus, in order to 
impose the Dirichlet boundary conditions we must con- 
sider the factor r a± . When Re(u>) — one can find that 



Re(a±) 



1 



± 2m, 



implying that one of ty± is convergent and the other one 
is divergent as r —¥ oo. The only exception is m = 0, 
when both and are convergent. 

For this case we recall that we require the function 
norm be convergent at spatial infinity, i. e., if R oc Ij^fr, 
then the norm 



J |$| 2 \/^d 4 xcx J \R\ 2 droz J y~m; 



(17) 



is divergent. The asymptotically anti-de Sitter space- 
times allow for a family of boundary conditions and a 
family of norms. For lower dimensional AdS black holes, 
such as the 2+1 dimensional BTZ black holes [l(|, the 
vanishing flux at infinity is used sometimes as the bound- 
ary condition [TlL fl2j . If instead of the above norm one 
uses the flux of the field, such norm will be even stronger 
divergent for the m = case. 

Thus, the boundary conditions do not allow us to have 
purely imaginary quasinormal mode for m = 0. However, 
we observe that Re(uj) approaches zero at some value of 
the rotation parameter (see Fig. [I}. Since we cannot find 
any quasinormal mode for larger values of the rotation 
parameter, probably, one should interpret this as disap- 
pearing of the particular mode after this threshold value 
of j is reached. 

We observe this "disappearing" of the quasinormal 
modes only for the case of m = 0. The quasinormal 
modes with the positive real frequency exist for all val- 
ues of j, if rn > (see Fig. [5]). For negative m we observe 
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FIG. 1: Real and imaginary parts of first three dominant quasinormal frequencies as functions of j (M = 1, fj, = 0) for A = 0. 
The real part of the dominant frequency is smaller and approaches zero for the largest threshold value of j. Higher overtones 
have larger oscillation frequency and damping rate. 

Re(w) 




FIG. 2: Fundamental modes as functions of j for I = m = 1 (blue, bottom), £ = m = 2 (green, middle), £ = m — 3 (red, top) 
(M = l, M = 0). 



Re(w) 




FIG. 3: Fundamental modes as functions of j for £ = 1 m = —1 (blue, bottom), m — (green, top on the left pane), m = 1 
(red, top on the right pane) (M = 1, /x = 0). 
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FIG. 4: Fundamental modes as functions of j for £ — —m — 1 (blue, bottom), I — — m — 2 (green, middle), I = — m = 3 (red, 
top) (M = 1, (j, = 0). 



the same tendency to "disappear" , though for the funda- 
mental mode the extremal value of j is reached before 
the real part vanishes (see Figs. G2 HJ). Since for higher 
overtones, modes disappear at lower threshold values of 
j, we may expect that for some near-extremal value of 
the rotation parameter we have only finite number of the 
modes with positive real part. Because of the symmetry 
of the spectrum, 

rn -m, Re(u) -Re(u), (18) 

we still have infinite number of modes with negative value 
of the real part. 

V. NUMERICAL METHODS 

Here, we shall briefly relate two standard methods used 
for numerical search of quasinormal modes: the shooting 
method and Frobenius method. These two alternative 
methods were used in order to guarantee the validity of 
the obtained results. More detailed discussion of both 
methods can be found in [7j. 

Shooting method. Since at the horizon we require 
the purely ingoing wave (|15p . it is convenient to define a 
new function in such a way, 

/ r 2\ iur Q y/\ - 4/i(r )/2 
y(r) = ( 1 " 3/ x (19) 

that it becomes regular at the event horizon. 

We fix the wave- function norm so that y{ro) — 1. Sub- 
stituting fll~9l into (TTTj) and expanding the wave equa- 
tion near the horizon, we find that y'(ro). This gives us 
boundary condition at the horizon for any fixed w, which 
we use for the numerical integration of the equation (jTTJ) . 
At large distances we compare the result of our numerical 
integration with the large-distance asymptotic series ex- 
pansion (TT51) and find the coefficients C± by using the fit- 
ting procedure. Then, quasinormal modes can be found 



by minimizing C+{uj) or C_(w) depending on the sign of 
Re(ui). In order to check convergence of the procedure 
we check that obtained frequencies do not change within 
specified accuracy, if we increase precision of floating- 
point operations or the distance at which we use the fit. 

Frobenius method. This method allows us to find 
QNMs by solving numerically an equation with contin- 
ued fractions, which takes much less computer time for 
finding quasinormal modes. One can rewrite the wave- 
function as 

/ 2 \-iur Q v / l-4:h(r )/2 
R(r) = e^ r r Q T I 1 - -| I x z(r), 

(20) 

where z{r) must be regular at spatial infinity and the 
event horizon, once u is the quasinormal frequency. We 
choose "-" sign for Re(ujj) > 0, and "+" sign for 
Re(ujj) < 0. 

The function z(r) can be expanded into a series near 
the horizon 

*(r)=5>(l-S) ■ (21) 

After we substitute ([20|) and (|21[) into equation (fTT|) . we 
find the three-term recurrence relation for the coefficients 
en 

Co,i flj + ci t i <ii_i + ca,i = 0, i>l (22) 
00,101+01,100 = 0. (23) 

From the above, the coefficients of the three-term re- 
currence relation can be found in a closed form. Finally, 
we find the equation with the infinite continued fraction 
on the righthand side 

C2.ATCo,AT-l C2,Af-lCo,Af-2 C2,2Co,l 

ci,w+i ■ ■ • ■ = 

Cl,jV-l— Ci,JV_2— Cl,! 

C0,N+lC2,N+2 Co,jV+2C2,JV+3 
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which we solve numerically with respect to the quasinor- 
mal frequencies ujn [13]. 



VI. THE LIMIT OF PURE GODEL D = 5 
SPACE-TIME (WITHOUT A BLACK HOLE) 

Let us take the limit M — > in the wave equation (|lip . 
Then, k(r) = f(r) = 1, the tortoise coordinate r+ — r 
and the wave equation (jlip is reduced to a much simpler 
form 

^ + ((1 - 4J V)- 2 - 8jmu - - ,, 2 )) i? = 0. 

(25) 

Substituting A = 1{Z + 1) and 

_ f e-^ r2 r- 2£ - 1 / 2 p(r), fle(jw) > 0, 
- | e jW 2 r -2^-i/2 p ( r ^ i? e (jw) < 0, 

we obtain the Kummer's equation 

zp"(z) + (b-z)p'(z)=ap(z), (26) 

with respect to the new coordinate 

2jtur 2 , Re(juj) > 0, 
-2jujr 2 , Re(juj) < 0. 

Here b = — 21 and 

f -8j(i - m)cj - w 2 + m 2 , fle(jcj) > 0, . . 

The solution of the Kummer's equation, which produces 
regular at the origin (r — 0) wave-function ^(r), is 

p(z) oc z 1+2£ M(a + 2^+1,2^+ 2, z). 

Here M(A _B, z) is the generalized hypergeometric series, 
which has an irregular singularity at z = oo, unless A is 
a nonpositive integer. Thus, one has 

a = -N-2£, N = 1,2,3.... 

Finally, from (|27|) we find 



lu = ±4j \ N + £±m+ \ (N + £±m) 2 + 



(28) 

One can see an interesting features of the scalar waves 
in D = 5 Godel Universe: there is a minimal value of 
frequency for any multipole number £, which is u> = 8j 
for /i = (see Fig. EJ. Waves with lower frequencies can- 
not propagate in the Godel Universe. A similar property 
was observed by Hiscock for the four-dimensional Godel 
Universe in [Tij ]. 

There are both positive and negative frequencies in the 
spectrum. Negative ones were interpreted by Hiscock [141 ] 
as those propagating back in time, because of the allowed 



closed time-like curves. Indeed, let us associate a wave 
with a negative frequency as that propagating back in 
time. Then, due to the symmetry (|18j) . we can conclude 
that such a wave with a given frequency corresponds to 
the opposite value of the azimuthal number m, which 
is eigenvalue of the time-like Killing vector in the region 
with CTCs. This means that the wave, propagating back 
in time in the region near the black hole, also propagates 
back in time when it reaches the region with CTCs. 

We can see that the perturbation spectrum of the five- 
dimensional Godel space-time is qualitatively similar to 
the four-dimensional spectrum [14j, and resembles the 
spectrum of empty anti-de Sitter space-time [HI, [l6[ , as 
both consist of equally spaced normal modes. 



VII. SPECTRUM OF A BLACK HOLE IN THE 
GODEL UNIVERSE 

In general, quasinormal modes u> are complex and 
Re(ui) is the real oscillation frequency, while Im{uj) is 
proportional to the damping rate of a given mode. For 
normal modes of pure Godel space-times discussed above 
Im(uj) = 0, so that waves look like nondamping (stand- 
ing) waves in a confining box. Another essential prop- 
erty of quasinormal modes is their significance for sta- 
bility analysis: if all QN modes are damped, the black 
hole is believed to be stable, while existence of any grow- 
ing mode indicates gravitational instability. Analytically 
stability is guaranteed if the effective potential is posi- 
tive definite outside the black hole, which happens only 
in a few simplest cases, i. e. for the Schwarzschild back 
hole. More general cases, and the Godel black hole is one 
of them, require further numerical approach to the wave 
equation for the analysis of of stability. Quasinormal 
modes have been effectively used in a number of works 
for checking of the gravitational stability (or finding of 
instability) of various black holes |17j . 

The pure Godel space-time has one scale parameter 
j. This parameter defines a scale of the universe which 
is distance oc j _1 to the region where the closed time- 
like geodesic trajectories exist. When we add a black 
hole to the Godel Universe, another scale parameter ap- 
pears, which is the black hole radius Tq. In this paper we 
consider fixed black hole masses and various universe's 
rotation parameters j. However, instead one could fix 
the scale of the universe j and look at various ro- In the 
latter case, black holes are parameterized by the dimen- 
sionless quantity Mj 2 . For small values of j, Mj 2 quickly 
approaches zero and the metric (Q| tends to the pure 
Godel space-time with the quasinormal spectrum given 
by the normal modes of Godel space-time (|28|) . what is 
confirmed by modes shown on Fig. [5j 

In concordance with this, on Fig. [S]it is shown that as 
M -> 0, that is Mj 2 -> at a fixed j, QNMs approach 
normal modes of the Godel space-time. This behavior ex- 
actly resembles the one of the asymptotically AdS black 
holes, for which QNMs approach normal modes of AdS 
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FIG. 5: First four dominant quasinormal frequencies as functions of j (M = 1, fj, = 0) for A = 0. As j approaches zero 
Re(uj) ~ 8Nj and Im{ui) oc j 4 . Higher overtones have larger oscillation frequency and damping rate. 

ReM 




FIG. 6: Fundamental mode approaches 8j as M —¥ (j = 0.1, X = 0, fi = 0). 



space-time in the limit ro — > fl5| . As j approaches zero 
at a fixed M, we have the limit of normal modes of the 
Godel space-time 

Re(uj) ~ 8iVj, Im(u>) oc j 4 , j -> 0. (29) 

Thus, at small j QNMs of black holes in the Godel 
Universe are shorter lived and have larger real oscillation 
frequency, when j increases. 

The numerical search of modes includes a kind of inde- 
terminacy because a mode is a root of quite complicated 
and lengthy algebraic equation, so that ad hoc we do not 
know which mode corresponds to which overtone number 
N. At the same time, we are most interested in funda- 
mental modes N = 1, which are longest lived and are, 
therefore, dominate at late time. In order to track value 
of TV of each mode, we started from very small values of j 
for which the quasinormal modes are close to the normal 
modes, and increased gradually the rotation parameter 
j- 



The Universe naturally contains only very small rela- 
tively the Universe scale black holes. Therefore, quasi- 
normal modes of asymptotically Godel black holes are in- 
evitably very close to the normal modes of the Godel Uni- 
verse. In other words, the quasinormal spectrum which 
we obtained represents rather the ringing of the Universe 
itself slightly corrected by the presence of a black hole. A 
somewhat more phenomenological approach to the black 
hole in the rotating Universe was suggested in @, in 
the limit of small universe's rotation. In Q the asymp- 
totic of the wave equation is qualitatively the same as in 
the 5-dimensional Schwarzschild case and the spectrum 
is Schwarzschild-like, though with some corrections pro- 
portional to j. 

Another interesting phenomena can be noticed on Fig. 
[I] Re(ui) vanishes after reaching some critical, moderate 
value of j, so that the mode cannot be numerically de- 
tected at larger j. As we could not fully investigate the 
limit of quasiextremal rotation j numerically, below we 
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shall try to analyze the wave equation analytically in this 
limit. 

QNMs which are shown on Fig. 1-6, although ob- 
tained for small and moderate j, do not demonstrate 
any tendency to instability, because the damping rate, 
which is proportional to Im(ui), increases as j grows. 
We have not reached, however, the extremal values of j 
for all modes and, therefore, cannot state that there is 
no instability for all j. 

VIII. THE LIMIT OF THE QUASIEXTREMAL 
BLACK HOLE IN THE GODEL UNIVERSE 

In the near-extremal black hole limit one has M — > 
l/8j 2 , while the event horizon approaches zero. Imply- 
ing that this situation is rather exotic, let us consider a 
universe rotating with speed very close to its extremal 
one, so that the radius of a black hole will still be finite 
and very close to zero. In this case, k(r) — > 1 and the 
tortoise coordinate r* again coincides with r. Then, the 
wave equation takes the following form 

d 2 R ( m 2 , 2 2 2 o i + 4A 2 \ „ „ 

(30) 

When m 7^ equation (|30[) has two irregular singular 
points: r — and r — oo, but when to = 0, r = 
becomes a regular singular point and the general exact 
solution can be found. For the massless field (// = 0) it 
is given by 

R = C\V^I- e -i/2^jr 2 ) + C 2 ^ + i/2^r 2 ), (31) 

where I a (z) is the modified Bessel function of the first 
kind. One of these Bessel functions converges at infin- 
ity, satisfying the Dirichlet boundary condition, while the 
other one diverges. However, this function is divergent 
also at r = 0. It is natural to expect that the function 
is divergent at the horizon. The radius of the horizon 
approaches zero when considering quasiextremal black 
holes. In order to have physically meaningful solution 
one has to require purely ingoing wave at the horizon, 
which is not possible for the extremal limit. Thus, for- 
mally the problem of the adequate boundary conditions 
in the extremal case remains open. If we require the 
wave-function to have a finite norm, there are no quasi- 
normal modes allowing it, what is consistent with the 
numerical results in the near-extremal region. 

Since (|3T)|) coincides with (|25l) for to = and fi 2 — > 
/i 2 — id 2 we repeat the procedure of Sec. I VII and find that 

8^a=(-^ + ^ W. U \>°> (32) 
\ Sjiu-^, Re(juj) < 0, v ' 

what cannot be satisfied for /i 2 > because a + £ < 0. 

For the particular case £ — 0, fi = 0, we can easily 
find the solution which satisfies the Dirichlet boundary 
conditions at spatial infinity 

R(r) oc e ±jur2 /y/r. 



TABLE I: Similarities in QN spectrum of asymptotically AdS 
and Godel space-times. 



asymptotically AdS BHs 
AdS radius R 


asymptotically Godel BHs 
Universe's scale j -1 


small BH (r < R) 
~ normal modes of AdS 


small BH (r < j -1 ) 
~ normal modes of Godel 


QNMs are equidistant at large N 



However, similarly to (|17p. one can prove that this solu- 
tion leads to a divergent norm as r — > 0. 

Thus, we conclude that for to = there is no so- 
lution, which satisfies the quasinormal boundary condi- 
tions. Numerical calculations show that this takes place 
in some near-extremal region as well. Physically we may 
interpret this as nonexistence of the most symmetrical 
perturbations of the quasiextremely rotating Godel Uni- 
verse, which can still be considered small. In case of 
arbitrary values of m and £ the question of stability of 
perturbations remains open. 



IX. DISCUSSIONS 

In the present work we have studied the scalar 
field perturbations around the Giomon-Hashimoto so- 
lution, which is the generalization of 5-dimensional 
Schwarzschild solution immersed in the Godel Universe. 
The main result of this paper is that the quasinormal 
spectrum of such "Schwarzschild-Godel" black hole is 
totally different from the Schwarzschild spectrum and 
resembles the spectrum of asymptotically anti-de Sitter 
black holes (see Table 13}. The latter is due to the Dirich- 
let boundary conditions at spatial infinity imposed for 
the asymptotically Godel space-times. The numerical 
analysis shows no signs of instability at small and mod- 
erate values of the Universe's scale j, which means that 
the conjectured correlation of the gravitational instabil- 
ity and existence of closed time-like curves, if it exists, is 
not straightforward. 

The frequencies U) are eigenvalues of the Killing vector 
dt which is timelike in some region near black hole. In 
the region far form the black hole we have CTCs and the 
timelike Killing vector is 9^ with the integer eigenvalue 
to. We cannot have a consistent solution which is not pe- 
riodic in the region with CTCs. Thus we can think about 
existence or nonexistence of growing oscillations only in 
the region near the black hole. The region with CTCs 
always implies periodical dynamics in it. This periodi- 
cal behavior of the scalar field in the asymptotical region 
effectively provides a "confining box" from the point of 
view of the observer near a black hole in the Godel Uni- 
verse. That is why we require the Dirichlet boundary 
condition at infinity. From the mathematical viewpoint 
this boundary condition is dictated by the structure of 
the irregular singularity at r = oo. 

In the future we would like to consider perturbations 
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of rotating analog of Gimon-Hashimoto solution, which 
should show much reacher physics. In particular, we ex- 
pect that the effective confining box given by the Dirich- 
let boundary condition together with rotation, and thus 
with superradiance, should produce the superradiant in- 
stability in the regime of quick rotation [18[ ■ 
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